Coarse-Grained Modeling of a Deformable Nematic Vesicle 



Jim Geng Jonathan V. Selinger, and Robin L. B. Selinger 

Liquid Crystal Institute, Kent State University, Kent, Ohio 44 24 2 > USA 
(Dated: December 18, 2011) 

We develop a coarse-grained particle-based model to simulate membranes with nematic liquid- 
crystal order. The coarse-grained particles form vesicles which, at low temperature, have orient a- 
tional order in the local tangent plane. As the strength of coupling between the nematic director 
and the vesicle curvature increases, the vesicles show a morphology transition from spherical to 
prolate and finally to a tube. We also observe the shape and defect arrangement around the tips of 
the prolate vesicle. 

PACS numbers: 82.70.Uv, 61.30.Jf, 87.16.A- 



Over the past twenty- five years, a major theme of re- 
search in condensed-matter physics has been the com- 
plex interaction of geometry with orientational order and 
topological defects. Both theoretical and experimental 
studies have investigated order and defects on surfaces 
of fixed shape, such as colloidal particles or droplets p]- 
[7]. These studies have shown, for example, that a ne- 
matic phase on a spherical surface will form four defects 
of topological charge +1/2 each, and these defects may 
be exploited to develop colloidal particles that will self- 
assemble into tetrahedral lattices for photonic applica- 
tions [8 j. Inspired by this potential application, many au- 
thors studied how to control the arrangement of the four 
half-charged defects in a nematic phase. Many effects 
have been considered, including elastic anisotropy [3], ex- 
ternal field [4 , and curvature of the colloidal particles [6] . 

Further research has investigated orientational order 
and defects on deformable vesicles, which serve as simple 
analogues for biological membranes [9HT2]. These stud- 
ies show that defects in the orientational order will de- 
form fluid vesicles into non-spherical shapes. For exam- 
ple, some vesicles have tilt order, which can be modeled 
as XY order in the local tangent plane; these vesicles 
will exhibit two defects of topological charge +1 each, 
and can deform into prolate or oblate shapes [11] . Other 
vesicles composed of T-shaped lipids or surfactants with 
rod-like heads may have nematic order in the local tan- 
gent plane [13]. Theoretical studies have predicted that 
the four half-charged defects will induce these vesicles to 
deform into tetrahedra [9]. 

So far, theoretical studies of orientational order in de- 
formable vesicles have considered systems that are ideal- 
ized in several ways: at zero temperature, with no elastic 
anisotropy, with only certain couplings between orien- 
tational order and curvature, and with shapes that are 
slight perturbations on perfect spheres. A key question is 
whether vesicle shape would be qualitatively different if 
any of these simplifying assumptions were relaxed. Com- 
puter simulation provides a useful approach to this is- 
sue. For example, simulations can investigate problems 
where the geometry is not a perfect sphere but rather 
a more complex disordered shape, with bumps of pos- 



itive and negative Gaussian curvature. One simulation 
method uses a triangulated-surface model with tangent- 
plane orientational order; this method has indeed shown 
complicated tube and inward tubulate shapes [12] . How- 
ever, a disadvantage of the triangulated-surface model is 
that the connectivity of the surface is fixed, unlike exper- 
imental systems in which molecules can detach from and 
rejoin the vesicles. 

In this article, we develop an alternative method to 
simulate orientational order in deformable nematic vesi- 
cles, using a coarse-grained particle-based model. This 
method allows the particles to be assembled into a mem- 
brane with orientational order in the local tangent plane. 
The membrane spontaneously selects its own shape, 
which may be flat, spherical, or more complex. Further- 
more, the interaction of the coarse-grained particles can 
be correlated with molecular features. Using this model, 
we calculate the arrangement of topological defects and 
the shape of the vesicle as a function of the interaction 
parameters. In particular, we find a morphology transi- 
tion from spherical to prolate and finally to a tube as the 
coupling between nematic order and curvature increases. 

To develop an appropriate simulation approach, we 
generalize a coarse-grained membrane simulation model 
without tangent-plane order due to Ju Li and collabo- 
rators P3H17; see also Ref. [15] . In their approach, a 
bilayer membrane is represented by a single layer of in- 
teracting coarse-grained point particles, each of which 
carries a polar vector degree of freedom n representing 
the preferred layer normal direction. The interaction po- 
tential favors association of particles with n vectors ly- 
ing parallel and side-by-side. In simulations, the particles 
self-assemble into pancake-shaped single-layer aggregates 
showing liquid-like self-diffusion and membrane elastic- 
ity. When another term is added to the potential favor- 
ing a slight splay between the n vectors of neighboring 
particles, particles spontaneously coalesce into spherical 
shells. Each particle represents not a single molecule but 
a large patch of membrane containing many molecules, 
and the surrounding solvent is implicit. This approach is 
more coarse-grained compared to other implicit-solvent 
lipid models in the literature [19j as it does not sepa- 
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FIG. 1. Each coarse-grained particle has a vector n, which 
aligns along the local membrane normal, and a vector c, which 
has nematic alignment within the local tangent plane. 

rately resolve hydrophilic head/hydrophobic tail compo- 
nents across the width of a lipid bilayer. 

To simulate a membrane with tangent-plane order, for 
instance formed by T-shaped lipids, we generalize the 
preceding model by defining a coarse-grained point par- 
ticle with two vector degrees of freedom, as shown in 
Fig. [I] The vector n again defines the preferred layer 
normal direction, and a new vector c represents the local 
nematic director orientation in the membrane's tangent 
plane. Both n and c are unit vectors and are always 
perpendicular to each other. Particles interact with each 
other via an anisotropic Lennard- Jones-type pairwise po- 
tential with a distance cut-off: 

^ij ( fi-i i i &i , C j , X^j ) = ( 1 ) 

u R (xij) + [1 + a[a(ni,nj,Ci,Cj,Xij) - i\]u A (xij). 

In this expression, the repulsive and attractive parts of 
the potential are given by 

ur{t) = e[(R cut - r)/(Rcut - r min )} 8 for r < R cut , 
ua{t) = -2e[(R cut - r)/(R cut - r min )] 4 for r < R cut , 
u R (r) = u A (r) = for r > R cut . (2) 

Note that the repulsive and attractive terms have expo- 
nents 8 and 4, respectively, in contrast with the expo- 
nents 12 and 6 for the Lennard- Jones potential. These 
reduced exponents soften the potential and enhance the 
fluidity of the membrane. The coefficient a controls 
the strength of the anisotropic orient ational interactions, 
which are defined by the function: 

^(tt-i) i Ci 5 Cj , X^j) (3) 

1 - [1 - (n, • n^ 2 - £]M(fi^) 2 " 7] 2 
-[(n,.x, J ) 2 - 7 ] 2 + 27 7 2 [(c i .c,) 2 -l]. 

In this function, /3 = sin 2 (6 ) o) and 7 = sin 2 (#o/2), where 
60 represents the preferred angle between the n vectors of 
two neighboring particles. As shown by Li and cowork- 
ers [14], if O =0? particles tend to self- assemble into flat 
membranes with the n vectors along the membrane nor- 
mal. If #0 7^ 0, then the n vectors of neighboring particles 
align with a separation angle #0; i- e - they favor a splay, 
and particles self-assemble into vesicles with spontaneous 
curvature. 



The novel aspect of our potential is the rj 2 term, which 
favors parallel or antiparallel alignment of the c vectors 
of neighboring particles, giving rise to nematic liquid- 
crystal order in the plane of the membrane. By varying 
the parameter n at fixed temperature, we can vary the 
strength of the nematic order, and hence the Frank elastic 
constants of the liquid crystal and the strength of the 
coupling between curvature and orient ational order. 

We perform a series of simulations with about 10,000 
particles over a range of n from 0.2 to 0.5, at tempera- 
ture ksT « 0.22e. Other parameters are d = 1.0, e = 1.0, 
a 3.1, r min 2 1 / 6 d, R cut 2.55d, and <9 = 0.0375 
radians. Initially, coarse-grained particles are distributed 
to cover the surface of a sphere of radius 30d using the 
random sequential absorption method [20]. The value of 
60 was selected to set the model membrane's spontaneous 
curvature to match the vesicle's initial radius. In the ini- 
tial state, n vectors are aligned radially and c vectors 
are randomly oriented in the local tangent plane. The 
temperature is increased and maintained by a Langevin 
thermostat. The c nematic director field quickly orders 
to form a population of +/— half-charge defects; this tex- 
ture coarsens via defect pair-annihilation until four posi- 
tive half-charge defects remain. The vesicle then relaxes 
toward its equilibrium shape. 

Varying the coupling parameter n gives rise to a se- 
quence of shape transitions. Snapshots of the front, side, 
and top views of the simulation results are shown in 
Fig. [2] Black dots indicate the locations of the four posi- 
tive half-charged defects in the tangent plane (except for 
the largest n = 0.5), and color indicates distance from 
the vesicle's center of mass. On the surface of the vesi- 
cles, the c vectors are shown as line segments in the local 
tangent plane. The n vectors are not shown. The coarse- 
grained particles are semi-transparent so that defects are 
visible on both the near and far sides of the vesicles. 

From the snapshots in Fig. [2] we can see that changing 
77 changes the nematic director field and defect configura- 
tion as well as vesicle shape. For the smallest n = 0.2, the 
four defects align themselves approximately on a great 
circle. This great circle is consistent with the contin- 
uum prediction of Ref. [3] for nematic vesicles with a 
strong elastic anisotropy, so it may indicate that this po- 
tential gives a substantial difference between the effec- 
tive Frank constants. As n increases to 0.25, the defects 
move to form approximately a regular tetrahedron. For 
these small values of 77, the vesicle shape is fairly close to 
spherical, indicating that the coupling between director 
distortions and curvature is not yet enough to substan- 
tially distort the sphere. As 77 increases further to 0.30, 
the defects shift further, with two defects moving to each 
end of the vesicle. At the same time, the vesicle elon- 
gates substantially between these two poles, forming a 
prolate shape. For n > 0.35, the nematic order becomes 
stronger, so that the defect cores cost more energy. We 
then observe holes nucleating at the centers of the defects, 



3 



r\ Front Side Top 




FIG. 2. (Color online) Simulated vesicles for several values of 
the parameter r\. Color indicates distance from the vesicle's 
center of mass, and black dots indicate particles in or near a 
defect core. Particles are semi-transparent so that defects on 
both near /far sides are visible. 



and the vesicle coexists with a particle gas. The defect 
configuration also changes from an elongated tetrahedron 
(rj =0.3-0.4) to rectangular (rj = 0.45). Finally, at the 
largest value 77 = 0.5, each pair of defects fuses to form 
a large pore at each end of the vesicle, thus eliminating 
the defects and transforming the vesicle into a tube. The 
middle of the tube is still swollen, not a perfect cylinder, 
because of the model membrane's spontaneous curvature. 




FIG. 3. (Color online) Enlarged front and top views of the 
vesicle for 77 = 0.3, with line segments showing the local ne- 
matic orientation c in relation to the overall shape. Particles 
are opaque; only defects on the near side are visible. 



To show the relationship between the local nematic 
orientational c and the overall shape more clearly, Fig. [3] 
shows front and top views of the vesicle with 77 = 0.3. In 
these views the coarse-grained particles are opaque, so 
that only two defects (on the near side) are visible, and 
the defects can easily be identified as topological charge 
1/2. From the front view, we see that the local nematic 
orientation c is aligned along the long axis of the prolate 
vesicle, allowing the c vectors of neighboring particles to 
be nearly parallel in three dimensions (3D); transverse 
alignment would have much higher 3D interaction en- 
ergy. From the top view, we see that the local nematic 
orientation c is aligned perpendicular to the separation 
between the two defects. Hence, the director distortion 
between the defects is almost entirely bend rather than 
splay. Furthermore, the top view of the vesicle is not 
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circular but extended along the average nematic direc- 
tor, perpendicular to the separation between the defects. 
Thus, the overall vesicle has a biaxial, potato-like shape. 

It is remarkable that the half-charged defects are ar- 
ranged in pairs, with a pair at each end of the vesicle, in 
spite of the usual repulsion between defects. We specu- 
late that this arrangement occurs because the defects are 
attracted to the regions of high positive Gaussian curva- 
ture at each end of the vesicle. This attraction to the 
curved regions competes with the repulsion between de- 
fects to favor an optimum separation between the defects, 
which depends on the coupling parameter rj. A similar 
pairing of defects has been seen in analytic calculations 
by Kralj [6] for the optimum positions of nematic defects 
on colloidal particles with a fixed ellipsoidal shape. Here 
we see that the pairing occurs even when the shape is not 
fixed but is free to deform. 

Our simulations can be compared with the predictions 
of Park et al |9] , who performed analytic calculations of 
the shapes of deformable membranes with nematic and 
general n-atic order (note n = 2 for a nematic phase). 
Their theory predicts that a nematic vesicle should de- 
form into a shape with the symmetry of a regular tetrahe- 
dron, with a half-charged defect at each vertex. By con- 
trast, our simulations never show the tetrahedral defor- 
mation, but only the extended biaxial, potato-like shape. 
We speculate that this discrepancy occurs because their 
theory is idealized in two ways. First, their free-energy 
functional considers only intrinsic coupling of director 
variations with curvature through a covariant derivative, 
explicitly neglecting other couplings allowed by symme- 
try [21]. The interaction potential in our simulation in- 
cludes an extrinsic coupling of the nematic director to the 
3D curvature direction. In previous work we showed that 
the extrinsic coupling greatly changes the director field 
on surfaces with fixed curvature [22]; here we see that it 
also affects shapes of deformable membranes. This effect 
may explain why the vesicle becomes extended along its 
long axis and why its ends become extended in a biax- 
ial way; both distortions reduce the interaction energy 
of c vectors in 3D. Second, their free-energy functional 
makes the approximation of a single Frank elastic con- 
stant, while our interparticle potential presumably gives 
different effective Frank constants. Shin et al [3] showed 
that anisotropy of Frank constants changes the arrange- 
ment of defects; here it also affects the membrane shape. 

The nucleation of pores in a bilayer membrane depends 
on processes at length scales well below the particle spac- 
ing in our coarse-grained model [23]. While our coarse- 
grained simulation results demonstrate qualitatively that 
pores spontaneously nucleate at topological defects, more 
detailed theoretical analysis and molecular-scale simula- 
tions are needed to understand the process in detail. An- 
other question is whether total enclosed volume is con- 
served during vesicle shape relaxation. The formation of 
pores in our simulation suggests that fluid may leak into 



or out of the vesicle as it evolves in shape, so we have not 
imposed a condition of fixed volume. 

In conclusion, we have developed a coarse-grained 
particle-based model for simulating membranes with ori- 
entational order, and we have used it to study vesicles 
in the nematic liquid-crystal phase. The simulation re- 
sults show surprisingly complex vesicle shapes and defect 
configurations, which arise from features in the interpar- 
ticle potential. Thus, the simulation method enables us 
to explore the range of phenomena that can occur in soft 
materials where geometry interacts with orientational or- 
der and topological defects. 

We thank A. Travesset for helpful discussions. Work 
was supported by NSF DMR-0605889 and 1106014. 
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